Shielding of absorbing objects in collisionless flowing plasma 

Yu. Tyshetski}tj and S. V. Vladimir ov 

School of Physics, University of Sydney 2006 NSW Australia 

(Dated: July 7, 2010) 

Abstract 

The electrostatic shielding of a charged absorbing object (dust grain) in a flowing collisionless 
plasma is investigated by using the linearized kinetic equation for plasma ions with a point-sink 
term accounting for ion absorption on the object. The efl'ect of absorption on the attractive part 
of the grain potential is investigated. For subthermal ion flows, the attractive part of the grain 
potential in the direction perpendicular to the ion flow can be significantly reduced or completely 
destroyed, depending on the absorption rate. For superthermal ion flows, however, the effect of 
absorption on the grain attraction in the direction perpendicular to the ion flow is shown to be 
exponentially weak. It is thus argued that, in the limit of superthermal ion flow, the effect of 
absorption on the grain shielding potential can be safely ignored for typical grain sizes relevant to 
complex plasmas. 
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I. INTRODUCTION 

In the field of complex plasmas [l-^, one of the central problems is the problem of 
plasma shielding of charged objects (dust grains) immersed in plasma. One of the effects 
that influences the shielding is the effect of absorption of plasma particles by a grain, with 
the latter acting as a sink for plasma particles, creating a flux of electrons and ions from bulk 
plasma, via the shielding region, to the grain. These fluxes affect plasma far from the grain, 
preventing it from reaching an equilibrium with the grain field, and hence affecting the grain 
shielding by the plasma in a qualitative way, compared to a non-absorbing (test) charged 
grain. This qualitative effect of absorption on the screening is important for interaction of 
dust grains in dust clouds, clusters and crystals , and for related phenomena such as phase 
transitions in complex plasma structures 6|. 

n 

It was recently shown [7| that in an isotropic plasma the shielding of a grain (that is 
at rest relative to the plasma) is changed qualitatively if absorption of ions by the grain 
is taken into account: instead of Debye-like potential at large distances from the grain, an 
tail appears in the potential, in case of coUisionless plasmas. However, in reality the 
plasmas in which dust grains are immersed are usually anisotropic, e.g., due to ion fluxes 
from bulk plasma to the plasma boundaries (these can be dust-plasma boundaries such as 
in dust voids [s], lol), or due to the motion of dust grains relative to the plasma (we note 
that the latter situation is equally relevant for satellites moving in the plasma environment 
of upper ionosphere 10|]). It is therefore necessary to investigate how absorption of plasma 
particles by a collecting body (e.g., a dust grain), either immersed in a flowing plasma, or 
moving in an isotropic plasma, affects the shielding of the body by the plasma. 



It is well known 
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12| that the potential of a test (non-absorbing) charge, either moving 



in isotropic plasma or, equivalently, stationary in a flowing plasma, has an attractive part 
in the direction perpendicular to the axis of the flow; in other words, two test charges of 
the same sign, placed sufficiently far from each other so that the line connecting them is 
perpendicular to the plasma flow direction, will attract each other electrostatically. This 
attraction may play a crucial role in formation and sustaining of plasma dust clusters and 
crystals (see and references therein). The attractive part appears due to the tail in 
the potential of the test charge far from the charge, which has the sign opposite to that of 
the test charge. It is reasonable to expect that, as in case of isotropic screening, absorption 
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will lead to a tail in the potential, of the sign opposite to that of the tail of the test 
charge potential, which will dominate over the tail at large distances from the grain, 
and hence may either reduce or completely destroy the attractive part of the potential, 
qualitatively changing the interaction of grains perpendicularly to the flow direction. 

The effect of absorption on grain shieldingjn plasmas with a flow has so far only been 
studied using a hydrodynamic model in Refs. 
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14j |. It was demonstrated that, at least in 



the case of highly coUisional plasmas considered there, the far asymptote of the grain field 
is dominated by the effect of absorption. However, it is yet not known how the absorption 
affects the grain shielding in plasmas with arbitrary degree of coUisionality, in particular, in 
weakly collisional or collisionless plasmas. To answer this question, a kinetic model of grain 
shielding, accounting for the absorption, should be used. However, the kinetic model for 
shielding of grains in anisotropic plasma developed by Ivlev et al. 15|], does not account for 



the absorption of ions by the grain. It is thus the aim of the present paper to study, using a 
kinetic model, the effect of absorption on grain shielding in the limiting case of collisionless 
flowing plasmas, thus covering the opposite limit of coUisionality to that studied in p^, Il4j |. 
In particular, we study here how the absorption affects the attractive part of the grain 
potential (in the direction perpendicular to the anisotropy axis) in collisionless plasmas. The 
answer to the latter question is important, e.g., for feasibility of the experiment proposed 
by Kompaneets et al. [l2| aimed at observing the attractive part of the grain potential in 
the direction perpendicular to the ion flow. 



II. MODEL AND GOVERNING EQUATIONS 

We consider a small spherical absorbing body (grain) of charge Q^, either moving in 
isotropic homogeneous collisionless plasma with velocity u, or immersed in a homogeneous 
collisionless plasma with a uniform flow u. Stationary kinetic equation for ions in the 
reference frame of the grain is 

^ ■ §^ - — V0(r) ■ ^ = -5{v)va{v)f, (1) 
or nii o\ 

where / = f{r,\) is the ion distribution function, (f){r) is the self-consistent potential of the 
grain in plasma. The term on the right-hand side (rhs) of ([1]) is the point-sink approximation 
of ion collection by the grain , in which 6{r) is the Dirac delta-function, and a{v) is the 
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ion collection cross-section of the body, which we assume to be isotropic, i.e., independent 
of the direction of ion velocity and only depending on its absolute value, o"(v) = cr(f). 
The assumption of isotropic ion collection cross-section is valid as long as the field of the 
grain, in which the ions being collected by the grain move, is approximately central, i.e., the 
anisotropy of the grain field in the flowing plasma is small in the region of size ~ ^ around 
the grain in which the trajectories of the ions being collected by the grain are significantly 
modified by its field. The rough criterion of validity of this approximation can be expressed 
as 



where is the characteristic length scale at which the anisotropy in shielding of the grain 
by plasma becomes significant. 

In collisionless plasma cr(f ) can be obtained from conservation of ion energy and angular 
momentum in the central field of the grain, and for small grains it is well approximated by 



the Orbital Motion Limited (OML) theory (l6| 



a{y) = aouM = Tra 1 , (3) 

where (pg is the surface potential of the grain, and a is the grain radius. Note that croML('^) 
does not depend on the distribution of potential in the viccinity of the grain. This is another 
reason for using the isotropic approximation a(v) = cr(f ) for ion absorption by the grain. 

Plasma electrons are assumed to be Boltzmann distributed, Ue = no exp(e0/Te), where 
no is the unperturbed plasma density (rieo = 'T-io = '^o); and is the electron temperature 
in energy units. The electron flux absorbed by the grain is assumed to compensate the 
absorbed ion flux, so that the grain charge Qd remains constant. The set is coupled by the 
Poisson's equation: 

- VV = 47re(ni-ne)+47rgd5(r), (4) 

where rii = J fdw is the ion number density, and (5d5(r) approximates the charge den- 
sity of the (small) grain, which is located at r = 0. This delta-function approximation is 
well justified when the grain size is small compared to the effective length scale of plasma 
screening. 

We employ the linear response formalism to find the static potential 0p(r) induced in the 
plasma by the absorbing charged grain. In the absence of the grain, the plasma is assumed 
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homogeneous and quasineutral, with no electric field (0o = 0), and the ion distribution 
function /o(v). The grain perturbs the plasma, inducing an electric field — V0p(r) in the 
plasma, and perturbing the distribution functions of ions, /(r,v) = /o(v) + /p(r, v). As- 
suming this perturbation to be small, \ fp\ <^ /o, we linearize Eqs solve them using 
Fourier transform, and obtain the self-consistent potential of the plasma perturbed by the 
grain in the form: 



Qd /■exp(zk-r) ie /" exp(ik ■ r) 

= 0QdW +0abs(r), 

where 0Qj(r) is the potential of the test (non-absorbing) charge Qd, and 0abs(i") is the 
additional term due to absorption of ions on the grain, in which 

*n...(k) = / l^^lMl),,. (6) 
J k ■ V — 

The dielectric function -D(k) in ([5]) is given by 

I,(,),l + %-4,/-5°W/f,v. (7) 



k"^ UQk"^ J k ■ V — zO 

where Upi = A/4vre%o7^ is the ion plasma frequency, and k^g = A^g is the inverse electron 
Debye length, A^e = VTe/^pe- The imaginary part — iO in ([7]) and ([6]) appears due to 
causality and specifies the rule for going around the pole when integrating over v. When 
the unperturbed ion distribution function /o(v) is approximated with a shifted Maxwellian 
in the reference frame of the grain. 



/o(v) = no$M(|v - u|) = ^^^^^exp [-^^ ] , (8) 
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where vxi is the thermal velocity of plasma ions, the dielectric function -D(k) becomes: 



. , TT /k ■ u\ ^ / k ■ u \ 



(9) 



2 \kVTi) ^ \^kvT^)\ ' 

where W{C,) = exp(— C^)erfc(— z^) is the plasma dispersion function of a real argument. 

Below we consider analytic solutions for 0p(r) of Eq. ([5]) in the limits of subthermal 
{u/vTi <^ 1) and superthermal {u/vxi ^ 1) fiow velocities. 
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III. RESULTS 



A. Subt hernial flows, u <^ vti 

This limit is relevant to a situation where a grain is suspended in a weakly anisotropic 
coUisionless plasma with a subthermal flow (e.g., a grain in coUisionless presheath), or to a 
grain slowly moving (compared to the plasma ion thermal velocity) in isotropic coUisionless 
plasma. 

Introducing C, = Is. ■ u/\/2kvTi and using the small argument expansion of W{C), we 
obtain for the dielectric function in the limit of subthermal flows: 

^(k) = 1 + § [1 - 2C' + ^V^C] + O(C^), (10) 

where fc^) = a/ + kj^e inverse total (electron and ion) Debye length. Note that in 

deriving ( |T0|) we assumed that electrons are hotter than ions, so that koe ^ koi and hence 
ko ~ koi- 

To calculate 5nabs(k) for u < Vxi, we use the approximation of the shifted Maxwellian ([8]) 
for small u/vti'- 

1 + ^ . (11) 

With this approximation, we have for 5nabs(k): 

/, N f vaiv) ^ , , , nnu f v\\va(v) ^ , , , 

y k ■ V — «0 f y k ■ V — 

Performing integration analytically where possible, we have (see Appendix |A]) 

5nabs(k) = eV(0 exp(-e^)rfe, (13) 

where ^ = v/\/2vTi- Finally, using ([3]) and integrating over ^, we obtain: 

5nabs(k) = ^7ra2(l + 2^r), (14) 

where r = Tg/Tj is the ratio of electron and ion temperatures, and z = e\(j)s\/Te is the 
dimensionless surface potential of the grain in plasma (here Te is the electron temperature 
in energy units). 

Inserting ( ITOj) and ( !T^ into ([5]), expanding up to the second order in u/vn, and employing 
spherical coordinates centered at the grain with the axis 9 = directed along u, we obtain 
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for and 0abs for u < vxi (see Appendix | 
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where 0^0 = eng ■ 7ra^(l + 22;r), and i?i and Ei are exponential integrals. Note that ( !T5|) 
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(16) 



matches the (somewhat corrected) solution of Cooper 11| for the potential of the slowly 



moving test charge in isotropic plasma. The structure of the field of the test (non- 
absorbing) particle has been analyzed in detail in [11]. Here we will only repeat some 
features of (pq^ relevant to this study, along with the corresponding features of the newly 
obtained 0abs- The spatial structure of the total grain potential 

Fig.m 

At small distances from the grain, we have for and 0abs: 



+ 0abs is shown in 



'abs 



exp(-fc£,r), 
r 

bao Mkor) +7-1] 



(17) 
(18) 



where 7 ~ 0.5772 is the Euler gamma. The ratio 0abs/0Qd tends to zero as r — )■ 0, i.e., the 
potential near the grain is dominated by the isotropic Debye potential f lT7|) . and the role of 
absorption on the near-grain potential is minor. 
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FIG. 1: Grain total potential (j)Q^ + (/>abs (normalized to Qd/^o), in case of subthermal flow, for 
different grain sizes (i.e., different rates of absorption): panel (a): k^a = (no absorption), panel 
(b): k£)a = 10~^, panel (c): koa = 2.5 • 10~^, panel (d): fczja = 10^^. The grain is located at 
p = z = 0, the ion flow u is directed from left to right. The = contour is shown with a thick 
black solid line on panels (a)-(c) [on panel (d) the potential does not change sign anywhere]. All 
potentials are plotted for u/vn = 0.2, r = 50, z = 3, no = 10^ cm~^. 



Outside the Debye sphere the behaviour of flTSj) is given by 



+ 



Q_ 

r 

2u 



vr u 



——klr^cos^e 
16 V 



kor > 1. 



(19) 



At large distances from the grain {knr ^ 1), the term proportional to in (fT9!) is dominant 
over the exponentially decaying term, and the asymptote of (pn at large distances from the 
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gram is 

, _ Qd 2m 
The far asymptote of 0abs of ( fT6ll is 



(20) 



•'abs 



Comparing the radial dependencies of and 0abs in ( 120]) and fl2T|) . we see that in case of 
subthermal flows, the absorption-induced potential 0abs dominates in the far-field potential 
of the grain, i.e., the far field of the grain in plasma is defined by absorption of plasma ions 
on the grain. [This result, obtained kinetically for collisionless plasmas, qualitatively agrees 
with the findings of Chaudhuri et al. 13(| for the dominant role of absorption in the far-field 
of grains immersed in strongly coUisional drifting plasmas.] The absorption of ions by the 
grain changes the radial dependance of the far asymptote of the total grain potential from 
oc r^^ (of a non-absorbing object) to oc r~^. The characteristic distance r(abs>Q) from 
the grain at which this change occurs is defined from ( l20l) and ( 1211) and depends on 9: 



\Qd\ u 2Xd u 

kDr{i,hs>Q} ^ ko— , 9 = 0,71, (22) 

(PaO VTi a VTi 

~i„M(^)'~H^(^)', e = ./2. (23) 

30aO \VTiJ 3a \VTiJ 

Here we have used z = e(f)s/Te with the surface potential of the grain approximated by 
0s ~ Qd/d, and the relation zt ^ 1 typical for dusty plasmas. 

The validity of the isotropic approximation for the ion absorption cross-section, as well 
as the validity of the linear theory expressions for the potentials and 0abs in case of 
subthermal flows, is discussed in Sec. IIVCI 



B. Superthermal flows, VTi u < Vs 



The limit VTi <^u < Vs (where Vs = {Te/niiY^'^ is the ion sound velocity) can be relevant, 
e.g., to dust grains suspended against gravity in a collisionless sheath region of the discharge 
(see, e.g., Q), or to satellites in the upper ionosphere 10|. In this limit, we can 

approximate the unperturbed ion distribution function with a shifted delta-function: 



/o(v) = lim no$M(v - u) = no5{v - u). 



(24) 




FIG. 2: Potential (normalized to Qd/^o) in the direction perpendicular to the ion flow, in case 
of subthermal flow. The thick dashed line shows the potential ()19p of a non-absorbing (test) 
charge, consisting of the Debye and r^'^ components shown with the thin solid and dashed lines, 
respectively. The thin dashed-dotted line shows the term (j2ip due to absorption of ions by the 



grain, and the thick solid line shows the total potential 



+ 



of the absorbing charged grain. 



The potentials are plotted for u/vxi = 0.3, a/A/) = 2 • 10 (panel a), and a/A^ = 10 (panel b). 



Neglecting the electron response (i.e., assuming that electrons are much hotter than ions so 
that electrons do not contribute significantly to the screening I4]), the dielectric function 
([9]) becomes in this limit: 



Z}(k) = 1 - 



1 



(A;|| -zO)2A2' 

where k\\ is the wavenumber along the direction of u, and A = u/upi = {u/vTijXm 
Integrating ([6]) with (!24|) . we obtain 



(25) 



UqU HQ 
C>nabs(k) = :j (J[u) = — (Toml(m) 

iC ■ U rh\\ 



(26) 



with (ToML defined in ([3]). 

Substituting (125|1 and (126|) into ([5]), we obtain for 



A 



00 



12| and for 0abs (see Appendix ICl) : 

z>0 (27) 



2i.o(0sin(| 



0Q.(P,-) = f / rft^^o (t^ exp (t^ 



z<0 (28) 
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FIG. 3: Plots of (l)Q^{p,z) and (j)ahs{p,z) given by Eqs ([27|) - (p0]) . Both potentials are normalized 
by the surface potential of the grain in vacuum, Qd/o,. The grain is located at p = z = 0, the flow 
velocity u is directed from left to right. The contours represent the equipotential lines (note that 
here the potential between the contours does not change linearly). The thick black contour shows 
the line of zero potential, separating the regions of positive and negative 4>/{Qd/a). 

and 

0abs(p,2;) = enoa{u) 



z>Q (29) 



where Jo is the zero-order cylindrical Bessel function of the first kind, Kq is the zero-order 
modified Bessel function of the second kind, and p, z are the cylindrical coordinates with 
the grain in the origin p = z = and axis z directed along u. Note that in writing out the 
Eqs ( 1271) and ( l29l) we omitted the common factor exp(— ■ z) that appears due to the term 
— iO in (125 p . which ensures that the field of the grain vanishes at 2; — -|-oo and finite p. The 
plots of (pQj{p, z) and 0abs(p, z) in case of superthermal flow are shown in Fig. |3l 

The potential fl27|) -f l28|) of the test charge Qd was obtained and analysed in detail in 12 |. 
At large distances from the grain, its asymptotic behavior corresponds to a quadrupole 
potential: 

^Q, = ^ (3 cos' ^ - 1) + 0(r-5), r ^ 00, ^ 0, (31) 
where 9 is the angle between vectors r and u. 
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The potential (129|) -( 130|) due to the effect of particle absorption on the grain has the 
following asymptotics: 

(jy^hs = enoa{u)^^^ + 0(r-^) cos^ 6, for 6 ^ 0, 7r/2, (32) 



-'abs 




enocr('u) W — exp I — — ) , for = 7r/2 (i.e., 2 = 0). (33) 



As we see, taking into account particle absorption on the grain leads to a qualitative change 
in the far asymptotic of the grain potential = + 0abs: instead of the quadrupole 
potential ( 13T]) . it becomes the dipole potential ( l32l) . The distance rdip at which the dipole 
potential ( 15^ starts to dominate over the quadrupole potential ([21]) is 



A 



enQo{u)\ 



(34) 



cos^ ' '2 

Hence in the case of superthermal ion flow the absorption is qualitatively important and 
should be accounted for, unless rdip/A ^ 1. In the latter case, the qualitative change in 
the grain field asymptote occurs at distances at which the field itself is extremely small 
and is overwhelmed by other fields. Therefore, in case of superthermal flows, the effect of 
absorption can be safely ignored for grains for which rdip/A ^ 1. After some algebra on 
Eq (J34l) . the criterion for ignoring the absorption becomes 



1 + zrf|.j/u^ aA cos6' 
where z = e|0s|/Te is dimensionless grain charge, 0^ is the surface potential of the grain, a 
is the grain radius, r = T^/Ti is the ratio of electron and ion temperatures, A^j is the ion 
Debye length, and A = u/ojpi = {u/vTi)\Di- 

The validity of the approximations made and of the obtained potential in the superther- 
mal limit will be discussed in Sec. IIV CI 



IV. DISCUSSION 



A. Modification of attraction between charges of the same sign 

It is easy to see from the asymptotes (1201) and (13 ip of for sub- and superthermal ion 
flows, that at = 7r/2 the far-field potential has the sign opposite to the sign of Qd, which 
implies that a non-absorbing charged grain will attract other grains with charges of the same 
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sign, if they are placed sufficiently far away from the test grain in direction perpendicular to 
the flow u 12| . Let us analyze how absorption of ions on the grain modifies this attraction, 
in the limits of subthermal and superthermal ion flows. 



1. Subthermal ion flows, u <C vxi 

The corresponding attractive potential well of a test (non-absorbing) charge in case of 
subthermal ion flow is shown by thick dashed lines in Fig. HJ^a-b). If this potential well is deep 
enough to confine other dust particles, a crystal-like structure of grains or a grain cluster 
might form in the plane perpendicular to the ion flow. The distance from the grain to the 
bottom of the attractive potential well defines the characteristic inter-grain distance in such 
structure, and the depth of the well defines the "melting temperature" of such structure, 
i.e., the characteristic kinetic energy necessary for a captured grain to escape the well. [For 
example, for the potential well shown in Fig. [21 the "melting temperature" of the dust grains 
of the same size as the test grain is of the order ~ Q^a/XoY eV. For a/A/) ~ 0.1 we have 
Trf ~ 6 ■ 10~^ eV, i.e., a dust grain with characteristic kinetic energy below 0.06 eV will be 
trapped in the attractive well of radius a ~ O.IA^ in the plane perpendicular to the ion flow, 
at distance from the test grain.] 

However, as can be seen from (1211) for 9 = it/ 2, absorption of ions by the grain introduces 
a term which decays as r~^, i.e., slower than the asymptote of oc —Qd/f^^ responsible for 
the attraction along the 9 = 71/2 direction, and has the sign opposite to it. In other words, 
depending on the parameters, absorption can significantly modify the attractive potential 
well of the grain, either making it much more shallow (and shifting its bottom further away 
from the grain) if r(abs>Q) > '^o, or destroying it completely if r(abs>Q) < '"o [here r(abs>Q) 
is defined by Eq. (I23p . and tq is the radial distance from the grain center to the bottom of 
its attractive potential well]. This modification is demonstrated in Fig. |21 It can be seen 
from Fig. [2] that absorption does not destroy the attractive potential well completely only 
for very small grains, while for larger grains absorption completely destroys the attractive 
potential well. Indeed, for u/vxi = 0.3 used in Fig. |2l the bottom of the potential well of 
the test charge Qd is at the distance Tq ~ 8A/) from the grain. The absorption will destroy 
the attractive potential well as long as the r"^ tail of 0abs starts to dominate over the 
tail of at distance less than the bottom of the potential well, i.e., if r(abs>Q) < ''^o- Using 
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(123|) . we find that this happens for grains larger than a/Xo ~ 10 ^, for the parameters of 
Fig. [21 For smaller flow velocities, u/vti < 0.3, absorption destroys the attractive potential 



well for even smaller grains. Typical range of dust grain sizes in dusty plasmas [19| is 
a ~ 0.3 — 20 yum, while typical range of Debye lengths is A/5 ~ A/^j ~ 20 — 200 yum, hence 
the typical ratio a/Xo is in the range a/Xr) ~ 10^'^ — 1, which only marginally includes 
the extreme value a/Xm ~ 10~^ obtained above. Thus we can conclude that in the limit 
of subthermal ion flows, the attraction between non-absorbing (test) charged grains in the 
direction perpendicular to the flow is completely destroyed if absorption of ions by the grains 
is taken into account, except for very small grains a/Xo < 10~^. 

2. Superthermal ion flows, vxi <^ u < Vg 

In the limit of superthermal flows, the corresponding attractive potential well in the 
direction 6 = 7t/2 from a non-absorbing charged grain has the minimum located at r 2.2A, 
and its depth is ~ 0.039Qd/X [l^. Absorption of ions by the grain introduces the term (155]) 
in the grain potential asymptote. This absorption-induced term is exponentially small and 
thus never dominates over the quadrupole asymptote fl3Tl) of 0q^. Thus we can conclude 
that in case of superthermal ion flows the effect of absorption does not destroy the attractive 
part of the grain potential in the direction perpendicular to the ion flow, for any grain 

his important conclusion is favorable for the experiment 



size and ion absorption rate. ^. 
proposed by Kompaneets et al [12| on detecting the attractive part of the grain potential 
in direction perpendicular to the flow. This experiment is to be performed with heavy dust 
particles levitated deep in the coUisionless sheath, for which the limit of superthermal ion 
flow, considered here, is applicable. 

B. Significance of the effect of absorption on the grain shielding potential 

Let us now outline the ranges of parameters for which the effect of absorption on the 
grain shielding is signiflcant and should be accounted for, and for which it can be safely 
ignored, in the limits of subthermal and superthermal ion flows. 
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1. Suhthermal ion flows 



Fig. [T] illustrates how absorption modifies the grain potential: as the grain size and the 
corresponding absorption rate increase, the attractive region around the grain shrinks, nar- 
rowing down to a finite-size region behind the grain, and eventually disappearing completely 
for large enough grains. The characteristic distance r(abs>Q), at which the absorption- induced 
term 0abs starts to dominate over the asymptote of the test charge potential 0q^, de- 
pends on and is given by Eqs (l22l)- (!23l) . The distance r(abs>Q) is the smallest along the 
6 = 7r/2 direction. If the minimum distance r(abs>Q) is not too large, r(abs>Q)(^ = 7r/2) < Tq, 
where is the distance to the bottom of the attractive potential well in the direction 
6 = 7r/2, then, as discussed in subsection IIV A T\ and seen from Fig. [H the effect of absorp- 
tion is significant for all grain sizes in the typical range a/Xo ~ 10^'^ — 1: for grain sizes in 
the range a/Xa > 10^^ the absorption destroys the attractive part of the potential along the 
9 = tt/2 direction, and for smaller grains the absorption modifies the attractive potential 
well along the Q = -n jl direction, making it more shallow. Hence, for grains of all sizes 



relevant to complex plasmas 



19| (a/A/) ~ 10 — 1) immersed in a collisionless plasma with 



subthermal ion fiow, one should account for the absorption-induced correction 0abs to the 
grain shielding potential, if one is interested in effects associated with the far asymptote of 
the grain field in the direction perpendicular to the fiow (e.g., grain interactions and phase 
transitions in dust clouds, clusters, and crystals). [We should note, however, that the ap- 
proximation of the grain as a point charge and point absorber in Eqs ([T]) and (jl]) becomes 
formally invalid for large grains, a/A/) > 1, in the subthermal limit. In this case the point 
approximation can be trivially generalized to account for the finite size of the grain. How- 
ever, such modification of the model will not affect the far field of the grain, and hence is 
not considered here.l 



2. Superthermal ion flows 

In the limit of superthermal ion fiows, one can safely neglect the effect of absorption on 



the grain potential if the criterion (!35|) is satisfied. Typical parameters [19| relevant to the 
case of superthermal ion fiows considered here (e.g., experiments in which dust grains are 
suspended against gravity in plasma sheath region with superthermal ion fiow) are: 2; ~ 1 
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(see Appendix [D]), r ~ 10 — 100, and u/vxi ~ ub/vti ~ 3 — 10 (here ub is the Bohm 
velocity of ions in the sheath, which is equal to the ion sound velocity Vg = ^/Te/rrii, and 
thus u/vTi ~ Vs/vti = \fTejTi = y/r). The criterion (|35i) then becomes, roughly: 

!iiiE ^ » 1, (36) 

A a 

which is well satisfied for typical grain sizes a/X^i ~ 10~^ — 1. 

We therefore conclude that for typical conditions of Earth-based dusty plasma experi- 
ments in which dust grains are suspended in the sheath region of the discharge with su- 
perthermal ion fiow, it is safe to neglect the effect of absorption of ions on the grain when 
calculating the potential of the grain in plasma, and just to use the expressions f l27|) - fl28|) 
for the potential of the test (non-absorbing) charge. Only for very large grains, a/Xoi ^ 10, 
suspended in a sheath plasma with superthermal ion fiow, one should account for the effect 
of absorption on the grain potential, given by Eqs (I29l)-(l30l). Yet we would like to stress once 
again that, as discussed above in Sec. IIV A2| the attractive part of the grain potential in 
the direction perpendicular to the superthermal ion fiow is only slightly affected by the ab- 
sorption, and can not be destroyed by it, for any grain size (unlike in the case of subthermal 
fiows) . 

C. Validity of the model 

Finally, let us discuss the validity of the approximations made in deriving the grain 
shielding potential in both subthermal and superthermal limits, and define the range of 
parameters for which the obtained results are valid. 

1. Linear approximation 

The total grain potential is calculated using the linear perturbation analysis, therefore 
the validity of the obtained expressions for the potential is limited to the region where the 
nonlinear effects are negligible, i.e., to distances sufficiently far away from the grain, such 
that the absolute value of the potential energy of ions in the grain field |e0(r)| is small 
compared to their characteristic kinetic energy e^. In the region near the grain where the 
potential energy of ions in the field of the grain becomes comparable to their characteristic 
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kinetic energy, |e0(r)| > ej, the shielding becomes nonhnear |20l-|22j . and this region of 
nonhnear shielding can be considered as a large "grain" of charge Qes {Qcs < Qd because of 
partial grain shielding in the nonlinear shielding region), which is then shielded linearly by 
the surrounding plasma. In this case the expressions for the grain potential still remain valid 
outside the region of nonlinear shielding, but with the true grain charge Qd replaced with the 
effective charge Qcs- However, if the characteristic size of the region of nonlinear shielding 
is small compared to the characteristic shielding length, then its effect can be neglected, and 
Qcs ~ Qd- 

Let us estimate the characteristic sizes of the nonlinear screening regions around the 
grain, corresponding to the two independent parts of the total potential of the grain, 
and 0abs- At small distances from the grain, the potentials (pg^ and 0abs in both subthermal 
and superthermal limits become 



Qd 

r 

4>aO 



ln(^)+7 



r < A„ (37) 
r < A„ (38) 



where 0ao = enocr(min{M, f^j}), A^ is the corresponding screening length [in the subthermal 
limit we have Xs = Xd and (pao = ^i^o<^{vTi), in the superthermal limit we have A^ = 
A = u/upi and (pao = enoO"(u)], 7 ~ 0.5772 is the Euler gamma. The characteristic sizes 
of the nonlinear screening regions around the grain, corresponding to and 0abs; are, 
respectively, r^^ ~ eQd/ei and r^^ ~ A^ exp[l— 7 — ei/(e0ao)] (note that the given expression 
for r^^ is vahd for r^^/As < exp(l — 7) ~ 1.53). 

In the subthermal limit, = Tj and A^ = Xd, and using fl37j) - fl38l) . we have for r'^^ and 



ZT, (39) 
AD Ad 

^ ~ expfl-7 (40) 

For the typical parameters 2; ~ 1, r ~ 10^, the distance r^^ is small compared to Xd only 
for very small grains, q/Xd < 10~^, while for larger grains the region of nonlinear screening 
is comparable to or even exceeds Ad, i.e., the shielding of the test (non-absorbing) charge 
is linear only for very small grains. However, as noted above, even for larger grains, for 
which r^^ ^ Xd and the shielding is nonlinear for r < r'^Ly linear theory still becomes 
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applicable at r > r^^^, and the expression f|T5l) for the shielding potential of the test charge 
is still valid for r > r^i, with Qd replaced by the effective charge Qes < Qd-, which is the 
remaining unscreened charge of the grain at the distance r^^ from the grain center. At the 
same time r^j^ is small compared to Xo for all grains with a/A/) < 0.1, while for a/A/) ~ 1 
(the largest grains of the typical range of a/A/) ~ 10^'^ — 1) rj^^ becomes comparable with 
Ad, and we expect some nonlinear modification of 0abs- Hence, the expression (1161) for the 
linear potential 0abs is valid for grains with a/A/) < 0.1, in subthermal limit. Yet even for 
large grains a/ A/) ~ 1, when the linear theory ceases to be valid, we do not expect the 
attractive part of the grain potential, already destroyed by the effect of absorption for much 
smaller grains a/A/) < 10^^, to reappear due to the nonlinear modification of (pahs- 
In the superthermal limit, using we have for r^^ and r^^: 

-S. '"^'--^Tl,r, (41) 



A Ae,- A e 



!:k_.p(i-.-?^(i..|)-). (42) 

For ej = 2 eV (in the sheath, to which the superthermal limit is relevant, ions can have 
even larger energies), Tj = 0.01 eV, Tg ~ (in the sheath, -u ~ t>s = ^/Tjrrii, and 
hence Te/e^ ~ 1), uq = 10^ cm~^, ~ 10^, and ^ ~ 1 (see Appendix [D]), we have 
r^^/A ~ 2;a/A ~ {a/ X£)){vTi/u) <^ 1 for a/A/) ~ 10^^ — 1, i.e., the linear theory expressions 
fl7r|) - fl25]) for (pQ^ are valid for all the grain sizes typical for complex plasmas [l9|. As for the 
ratio r'^il A, it is extremely small for the wide range of parameters, including those typical 
for the sheath region where the superthermal limit is applicable, due to the exponent in 
fH2|) . Hence, the linear expressions fl29|) - fl30|l for 0abs in the superthermal limit are valid for 
all of the practical grain sizes. 



2. Isotropic approximation for ion absorption 

In our model, we also made the approximation of isotropic cross-section of ion absorption 
on the grain, cr(v) = cr(f ). The validity of this approximation is defined by the criterion (jj]). 

In case of subthermal ion flows, the distance A^ at which the anisotropy of the grain field 
becomes significant is roughly defined as the distance at which the part of (|T9l) starts 
to dominate over the Debye part of (fT9l) . which is of the order of several A/) for u/vxi < 1, 
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and increases with decreasing u/vn (for u/vn = 0.1 we have Xa ~ 5.8Ad, and for u/vn 
the shielding becomes isotropic, and Xa — ?■ oo). Taking a ~ 7ra^(l + 2zr), we obtain from 
([2]) the following criterion for validity of the isotropic approximation: 

Xa/Xo 



< 



Xd VI + 2zt' 



(43) 



For z ~ 1, r ~ 10^, and ~ lOA/j, the isotropic approximation cr(v) = cr(f) is valid for 



a 

Ad 



(44) 



i.e., for most of the typical grain sizes o/Ad ~ 10^'^ — 1. Even for the largest of the typical 
grains, since in the OML approximation a{v) does not depend on the distribution of the 
grain potential in the vicinity of the grain, we do not expect a significant effect of the grain 
field anisotropy on the process of ion absorption on the grain, and hence on the derived 
absorption-induced term 0abs in the grain potential. 

In case of superthermal ion flows, we use ([2]) with A^ > A, cr = cr{u), which gives us the 
following criterion of validity of the isotropic approximation for ion absorption on the grain: 

u/vTi 



< 



(45) 



Ad ^yl + zTe/ei 

For z 1, u/vTi ~ y/r ^ 1, and Tg/ej ~ 1, this condition is well satisfied (and hence the 
isotropic approximation is valid) for all grain sizes from the typical range a/Ad ~ a/Xoi ~ 



10-3 - 1 
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We therefore conclude that the approximation of isotropic cross-section of ion absorption 
on the grain is well justified for typical grain sizes relevant to complex plasmas. 

D. Effect of ion-neutral collisions 



The model used in this paper is collisionless, but it can in principle be generalized to 
include ion-neutral collisions, which are typically the most frequent collisions in complex 
plasmas [16|. The effect of ion- neutral collisions will most probably modify the results 
obtained here, as the collisions lead to an tail in the shielding potential of a test (non- 
absorbing) grain in drifting plasmas 23|, l24|^ and to an tail in the shielding potential 



of an absorbing grain in isotropic plasmas 



7 
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26] or drifting strongly coUisional plas- 



mas [13j. The generalization of the kinetic model with absorption used in this work, to 
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include collisions, and the study of the effect of absorption on shielding of dust grains in 
flowing plasmas with arbitrary degree of ion-neutral coUisionality, are left for future work. 



V. CONCLUSIONS 

In this work we obtained the shielding potential of an absorbing charged grain, which is 
either moving in a uniform coUisionless plasma, or is immersed in a collisionless plasma with 
a uniform flow, employing linear perturbation analysis of the kinetic equation for ions with 
point-sink term accounting for the absorption of ions by the grain. The total potential of the 
grain in plasma is shown to consist of two independent terms: the term (pq^ due to the grain 
charge, and the term 0abs due to absorption of plasma ions by the grain. Having obtained 
the general expression for the absorption- induced term 0abs, we analyzed analytically its 
effect on the grain shielding in the limits of sub- and superthermal ion flows. It has been 
shown that, in general, absorption leads to a qualitative change in the far field of the grain 
from a quadupole-like field of a non-absorbing (test) charge to a dipole-like field of the 
absorbing grain. The characteristic distance from the grain at which this change occurs 
is obtained in the limits of subthermal and superthermal ion flow in terms of the plasma 
parameters and the grain size. This distance varies depending on the direction relative to the 
direction of the flow, and so does the effect of absorption-induced correction 0abs to the grain 
shielding potential. The most important absorption-induced effect on the grain shielding is 
the effect on the attractive part of the grain potential in the direction perpendicular to the 



flow, t 
grain 



lat has been previously demonstrated to appear for a test (non-absorbing) charged 
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12|, and which is important for the grain interaction in structures (e.g., dust 
clusters and crystals) in complex plasmas. 

In the limit of subthermal flows, the absorption of ions by the grain is shown to signif- 
icantly modify this attractive part of the grain potential. The character and magnitude of 
this modification depends on the plasma parameters and the grain size. The quantitative 
criterion is derived which defines, for given plasma parameters and grain sizes, how the ab- 
sorption modifies the attractive part of the grain potential: whether it makes the attractive 
potential well more shallow, thus making the attraction weaker and reducing the associated 
"melting temperature" of the dust structure in the plane orthogonal to the flow direction, 
or whether it destroys the attractive part of the grain completely. Our estimates show that 
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in the subthermal limit u/vxi <^ 1 the effect of absorption on the attractive part of the 
grain potential is significant for all of the typical grain sizes a/Xo ~ 10^'^ — 1 relevant to 
complex plasmas. In particular, we found that in case of subthermal flow the absorption 
significantly reduces the attractive part of the grain potential for the smallest of the typical 
grains, a/Xo ~ 10^^, while completely destroying it for grains larger than a/Xn ~ 10~^. 
One should thus account for the effect of ion absorption on the grain in plasmas with sub- 
thermal ion flow, if one is interested in the effects associated with the asymptotic behavior 
of the grain potential in the direction perpendicular to the flow, e.g., interaction of grains 
and phase transitions in dust crystals and clusters in the plane orthogonal to the ion flow. 

In the limit of superthermal flows, however, the situation is different from that of the 
limit of subthermal flows: the absorption can only slightly modify the attractive part of the 
grain potential in the direction perpendicular to the ion flow, but can not destroy it, for 
any grain size. Moreover, for typical parameters and grain sizes relevant to complex plasma 
experiments where dust grains are levitated in the sheath region, the minimum distance at 
which the absorption-induced change from the quadrupole to the dipole field occurs is so 
large compared to the characteristic screening length A [see Eq. ( I36l) ]. that the grain field at 
this distance is practically shielded completely, or at least is shielded to the level where it is 
overwhelmed by random plasma noise fields. Therefore, in case of superthermal ion flows, 
the effect of absorption on the grain shielding potential can be safely ignored for typical 
grain sizes a/Xoi ~ 10"^ — 1 occuring in complex plasmas. 
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Appendix A: Derivation of Eq (I13p 

Using d\ = v±dv±dv\\da, the integral in the second term in the right-hand side of ( fT2|l 
can be written as 



/ 



v\\va{v) 



(^Miy)dw = 2lT 



oo 



(Al) 



oo 
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where 

G{vl)= r dv^v^^J^k^=<l>M{v±), (A2) 

with V = + v'j_. Since G is an even function of the above integral is zero, and the 
second term in the right-hand side of (fT2!) vanishes. 
The integral in the remaining term in flT2|) is 

= ,i„„ / ^^t„(„)<iv. (A3) 
J k ■ V - !^->o y z/ + ■ V 

Using spherical coordinates (iv = sm{9)dvd9dip and integrating over 6' and ip, we reduce 
the last integral to 

va{v) 



-^Miv)d\- = arctan exp(-^^)rf^, 



(A4) 



+ ik ■ V 

where ^ = v/ -\/2fTj, 6 = z// \/2kvTi- Taking the limit z/ — )■ in ( 1A3p . we obtain the Eq. (fT3|l . 
Appendix B: On derivation of Eqs (|15|) and (|16p 

To evaluate the integrals on k appearing after substituting f lTOj) and f lT4|) into ([5]), we 



expand exp(zk ■ r) following Cooper 



exp(ik-r) =47r5^ J2 t'ji{kr)Yil{ek,ipk)YU0,^), (Bl) 

/=0 |m|<« 

where ji are the spherical Bessel functions, Yim are spherical harmonics, 6, ip and ^fc, (pk are 
the spherical angular coordinates of the vectors r and k, respectively, and the asterix denotes 
complex conjugation. Then, substituting (IBip into (E]) and using dk = k"^ sin 9 kdkdOkdipk, 
we perform the integration over 6k and (pk and then, where possible, over k, obtaining the 
Eqs ([15]) and (dS]). 
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Appendix C: Derivation of Eqs ([271) - ([30l) 



Upon substituting (125|) and ( 126|) into ([5]), we obtain the following expressions for and 
0abs in cylindrical coordinates p, 2; (note that k ■ r = k^pcosipk + (ik = k±_dk^dk\\dipk)'- 

Qd^^ f°° f°° k?,exp{ik\\z) 
(t>Qd{p,z) = / dk±k±Jo{k±p) / dk\\- ^ , (CI) 

^ ^0 J-00 (kj + kij {{k\\ - ioyx^ - 1) 

ieriQ 2 / \ /"°° \ f°° k\\exp{ik\\z) 
4>abs{p,z) = A a{u) / dk±k±Jo{k±p) / dkny r ,(C2) 

^ Jo J-00 f A;2 + A;2 j ((A;|| - iO)2A2 - 1) 

where the integration over the azimuthal angle has been carried out using the formula 

pin 

I exp{ik±p cos (pk)d(pk = 27rJo{k±p). (C3) 
Jo 

The integration over k\\ in (ICip and (IC2p can be performed analytically using the residue 
theorem and the Jordan lemma. The poles of the integrand are k\\ = ±ik± and k\\ = 
±A~^ + iO, and the contour of integration is closed in the upper (lower) part of the complex 
plane of fcy for z > (2; < 0). Performing the integration over fcy and using the formula 

jr°°,,^%M^ = ,,„^^). (C4) 

we obtain Eqs (|271)-(15D|). 



A;2+A-2 ^VA 



Appendix D: Estimate of the grain charge in plasma with superthermal ion flow, 

U » VTi 

The grain charge reaches equilibrium when the electron and ion currents on the grain 
compensate each other, 1^ = h- In the superthermal limit, we approximate the ion distri- 
bution function with /i(v) = no5(v — u), thus we have for the ion current on the grain: 



va{v)fi{v)dv = riQuaiu). (Dl) 
Using the OML approximation ([3]) for a, we have 

li = n^uira^ ( 1 + 7-^^ ) • (D2) 



For Maxwell-Boltzmann distribution of electrons in the grain field, the electron current 



on the grain is 



16| 



Ie = V Sna noVTe exp ( — ) , (D3) 
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where 0s is the surface potential of the grain. 

Now, introducing the dimensionless quantity z = e|0s|/Te (normahzed grain surface po- 
tential), which is approximately equal to the normalized grain charge, z ^ e^Z^/aTg (here 
Zd is the charge number of the grain), and requiring that I^. = li, we obtain the equation 
for z: 

^exp{-z) = —.[^ (l + 2zT^] , (D4) 

VTi \ T \ J 

where /i = rrie/mi is the ratio of electron and ion masses, and r = Tg/Ti is the ratio of 
electron and ion temperatures. The solution z of flD4p weakly (logarithmically) depends on 
the plasma and flow parameters /i, r, and u/vTi, and is of the order of unity, z ~ 1, for 
typical values of these parameters. For example, for a grain suspended deep in the sheath 
region we have m ~ f^, where Vg = \fTeJrni is the ion sound velocity, and hence ujvTi ~ ^fr. 
Then the equation (]D4p reduces to 

exp(-2) ~ ./^(l + 2z). (D5) 
V 8 

For hydrogen plasma (;U 5.44 ■ 10~^) its solution yields z ^ 2.45, while for argon plasma 
(^ ^ 1.36 • 10~^) its solution yields z ^ 3.90. 
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